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Complete Genome Sequences of Five Paenibacillus larvae 
Bacteriophages 
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Paenibacillus larvae is a pathogen of honeybees that causes American foulbrood (AFB). We isolated bacteriophages from soil 
containing bee debris collected near beehives in Utah. We announce five high-quality complete genome sequences, which repre- 
sent the first completed genome sequences submitted to GenBank for any P. larvae bacteriophage. 
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Paenibacillus larvae is a facultative anaerobic spore-forming 
pathogen that causes American foulbrood (AFB). AFB kills 
honeybee larvae (1), contributes to colony collapse disorder (2), 
and limits agricultural yields (3). Unfortunately, some P. larvae 
strains have become resistant to the antibiotics typically used 
for AFB treatment (4). Phage therapy is a potential treatment 
for AFB, yet few P. Zarvae-specific phages have been described 
(5-7), and a full-genome sequence for one has only recently 
become available (8). 

Soil, honey, and larva samples were collected in the Utah, Salt 
Lake, and Davis counties of Utah. The samples were used for iso- 
lating new host strains or bacteriophages. P. larvae subsp. pulvifa- 
ciens (9) hosts were confirmed by 16S rRNA sequencing. The bac- 
terial cultures were inoculated with soil samples to enrich for 
bacteriophages. Plated enrichment samples formed plaques from 
which bacteriophages were selected and purified by a minimum of 
three passages. High-titer lysates were filtered, incubated with 
5 ju-g/ml RNase and 10 ;u,g/ml DNase for 30 min at 37°C, and 
treated with 100 /xg/jLtl proteinase K at 52°C for 1 h. Following 
phenol-chloroform extraction and ethanol precipitation, high- 
quality DNA (10) was sequenced using 454 pyrosequencing. 

Raw sequences were assembled into contigs using Newbler ver- 
sion 2.6 (Roche Diagnostics, Branford, CT) and Consed version 
19(11). The phages Abouo and Emery assembled into single con- 
tigs. The phages Jimmerl, Jimmer2, and Davies assembled into 
multiple contigs that were joined using Gepard 1.30 (12),MEGA5 
(13), and Geneious Pro 5.4.4 (Biomatters Ltd., Auckland, New 



Zealand) with phage Abouo as a reference. The sequencing fold- 
coverage data are provided in Table 1. 

No experiments were performed to determine the physical 
ends or the packing or replication strategies of the phage DNA. 
However, during manual finishing, overlapping contigs assem- 
bled the genome into an apparently circular genome. For the bac- 
teriophages Jimmerl, Jimmer2, Davies, and Abouo, the first base 
of each genome was selected in the noncoding gap between the 
terminase gene and the prior gene. Since Emery did not have a 
terminase small subunit gene, the first base was selected in the first 
gap upstream of the large subunit terminase gene. The noncoding 
gap where the first base was selected for a genome contained mul- 
tiple stops and lacked coding potential in any frame in all ge- 
nomes. 

Annotation was completed using DNA Master (http://cobamide2 
.bio.pitt.edu). A coding potential map was generated using 
GeneMark 2.5p ( 14) for each phage based on Bacillus cereus strain 
ATCC 14579, the closest available relative to P. larvae. Our 
selection of gene calls emphasized the following criteria: 
GeneMark HHM and Glimmer autoannotation, BLAST align- 
ment E values of <0.001, coding potential from GeneMark, start 
codon sequences, and Shine-Dalgarno (SD) scores of >200 nats 
using the Karlin position-specific scoring matrix (PSSM) for 
moderately to highly expressed genes. 

All five phages were identified as myoviruses. Jimmerl and 
Jimmer2 were isolated independently from the same soil sample 
and differ in their genomic sequences by only 80 bp of the 



TABLE 1 Paenbacillus larvae bacteriophage 


genomes 








Phage name 


GenBank accession no. 


Sequencing fold coverage 


Length (bp) 


No. of genes 


G+C content (%) 


Jimmerl 


KC595515 


250.2 


54,312 


100 


38.11 


Jimmer2 


KC595514 


271.5 


54,312 


100 


38.10 


Emery 


KC595516 


143.2 


58,572 


100 


41.44 


Abouo 


KC595517 


116.9 


45,552 


92 


39.16 


Davies 


KC595518 


130.8 


45,798 


93 


39.16 
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54,312-bp genomes (99.85% similarity). The differences between 
these related phages are real base pair changes because the se- 
quencing fold coverages of these samples are statistically greater 
than the error rate of 454 sequencing (15). 

Nucleotide sequence accession numbers. The GenBank ac- 
cession numbers for the five Paenibacillus larvae bacteriophages 
are listed in Table 1. 
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